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Abstract. The high-pressure dynamics of a computer-modeled silica melt is 
studied in the framework of the mode-coupling theory of the glass transition 
O ■ (MCT) using static-structure input from molecular-dynamics (MD) computer 

'^ ' simulation. The theory reproduces the experimentally known viscosity minimum 

(diffusivity maximum) as a function of density or pressure and explains it in 
terms of a corresponding minimum in its critical temperature. This minimum 
arises from a gradual change in the equilibrium static structure which shifts from 
being dominated by tetrahedral ordering to showing the cageing known from high- 
density liquids. The theory is in qualitative agreement with computer simulation 
C^ ' results. 



o 

,^,' 1. Introduction 



PACS numbers: 64.70.Pf, 62.50.-|-p, 66.10.-x 



The physical mechanisms behind vitrification are still widely debated. Studies taking 
^ ' into account the high-pressure behaviour of glass-forming liquids in addition to their 

^—^ . response to temperature variation, are now emerging as a valuable means to provide 

insight into the glass transition phenomenon. Based on a comparison of colloidal 

fvq . and molecular glass-transition data over large pressure ranges, it has been proposed 

that pressures well above 1 GPa are needed to significantly change the properties of 
many well-studied fragile glass formers [l|; only in this extreme pressure regime one 

QQ ' would have hope to resolve the long-standing debate whether energetic or entropic 

f^ , interactions are the main cause for the dynamical slowing down in the vicinity of the 

glass transition. 

It has been pointed out 0] that most organic glass formers will cease to exist 

/\ ' as such under these conditions, the molecules being irreversibly transformed to 

polymerized modifications. Even some of the existing measurements on organic glass 
formers, touching on the 1 GPa regime [Sj], might have to be reconsidered in light of 
this findingU This leaves two material classes as possible candidates for furthering 
the understanding of pressure-induced vitrification and, by this route, the physics of 
the glass transition itself: metallic glasses on the one hand, and amorphous silica and 
relatives on the other. 

I Brazhkin V V, private communication 



The Dynamics of Silica Melts under High Pressure 2 

Silica and silicate melts are special in that they are both ubiquitous in application 
and known to display 'anomalous' changes in thermodynamic and kinetic properties 
as pressure is increased. This arises essentially because they exist as open tetrahedral- 
network-forming structures under atmospheric conditions and can be classified as 
'strong' liquids in Angell's sense. A prominent feature in silica and a number of alkali 
silicates is that the diffusivity of the Si and O atoms first increases with increasing 
pressure \4, la, ly, lZ| , contrary to what one expects from the excluded- volume picture of 
the glass transition, where increased pressure, insofar as it leads to increasing density, 
drives dynamic arrest. Only at pressures higher than about 10-20 GPa, depending 
on the melt composition, do the diffusion coefficients decrease for silica melts, leading 
to a maximum in the diffusivity- versus-pressure plot, as predicted from computer 
simulation [Sl, M, llO| . Similarly, the viscosity first decreases with increasing pressure, 
eventually showing a minimum. 

The mode-coupling theory of the glass transition, MCT 'lli|, is commonly 
accepted to be a theory applicable to relatively high temperatures in 'fragile' liquids 
where it has been tested with great success [l^] , although MCT signatures such as two- 
step slow relaxation have also been seen in 'strong' liquids displaying nearly- Arrhenius 
behaviour at low temperatures [l3|. Earlier simulation and combined simulation- 
and-theory studies [ij, [iSj point out by demonstrating qualitative agreement with 
microscopic MCT calculations for silica and sodium silicate mixtures, that the theory 
can in fact yield more detailed predictions also for strong glass- forming liquids. The 
existence of a MCT-Tc also for silica implies a number of asymptotic predictions, for 
example that the long-time 'structural relaxation' in the system is in fact independent 
on the details of microscopic motion; this strong MCT prediction has recently been 
confirmed in a simulation study comparing molecular-dynamics with stochastic Monte- 
Carlo dynamics for a simulation model of silica |16| . Still, it remains a crucial question 
towards understanding the theory's benefits and limitations, how MCT deals with the 
differences between 'strong' and 'fragile' liquids. Silica has been argued to undergo a 
transition from one to the other upon pressurization [9] , and so is an ideal candidate 
for these studies. 

In this contribution, we present first MCT results for pressurized silica melts. 
We demonstrate that the diffusivity maximum is reproduced by the theory, using 
computer-simulated static structure factors as input. The maximum of diffusivity 
corresponds to a minimum of the MCT critical temperature Tc as a function of density 
g or pressure P, which in turn arises from the interplay of decreasing tetrahedral- 
network effects and increasing contributions from nearest-neighbour cageing. 

2. Model Calculations 

Since MCT predicts drastic changes in the dynamics arising from relatively minor 
changes in the average equilibrium structure, one needs to ensure good-quality input 
for the latter. For complex liquids such as silica, this can be delivered by computer 
simulation in conjunction with a reliable model potential. Carre et al |l7| have 
recently developed a pair potential (called CHIK potential) based on Car-Parrinello 
calculations that reproduces the experimental equation of state reliably and hence 
improves significantly over the Beest-Kramcr-van Santen (BKS) potential [l8( widely 
used so far. The CHIK potential was used in extensive molecular-dynamics (MD) 
computer simulations described in detail elsewhere in this issue |l9| . From these, 
the matrix of equilibrium partial static structure factors S{q) = {na{q,t)*nf3{q,0)) 



The Dynamics of Silica Melts under High Pressure 3 

was obtained, where na{q,t) — "^iGxpliq- ri^ait)] are the Fourier-transformed (wave 
vector q) number-density fluctuations of species a — Si, O, and the sum runs over 
all Na particles belonging to that type with positions fi^a{t). For use in the MCT 
equations, S{q) has been simulated for the isotherms T — 2100 K, 2230 K, 2580 K, 
2750 K, 3250 K, and 3580 K at various densities from 2.3g/cm^ to 4.3g/cm'^; at other 
temperatures, S{q) was obtained by linear interpolation between the above values 
unless otherwise noted. 

MCT equations of motion for the resulting binary mixture are solved numerically 
as outlined in Ref. [20|. They yield the matrix of dynamic partial number-density 
correlation functions $(q,i) depending on wave number q — \q\, 

J-\q)d!^{q,t) + S-\q)^{q,t) 

(1) 



/ M{q,t-t')dt'^{q,t')dt' ^0, 
Jo 



where Japiq) — q^kBT jmaSap sets the thermal velocities for the short-time dynamics 
and $((7, 0) — S{q). M{q, t) is the memory function matrix of generalized fluctuating 
forces, which in the MCT approximation is written as 

, ^ , N 1 n f d^k 

X ^ Vaa'a"iq,k)Vf3fj'P"{q,k)<^a'l3'{k,t)<^a"l3"iP,t) 
a'l3'a"0" 

with p = \q — k\. Here, n is the number density, and Xa are the number 
concentrations, a;si — 1/3 and xq — 2/3 in our case. The vertices VaQ'Q"(<f, fc) = 
{qk/q)caa'{k)Saa" + {qp/ q)caa" {p)Saa' + qnx ac)^^, ^„ {q, k) Contain only equilibrium 
static correlations, viz. the matrix of direct correlation functions c(q) defined by •S'(g) 
through the Ornstein-Zernike relation [2l| . c^^^ denotes the corresponding static triplet 
correlation function [22| . which we set to zero in the following since no simulation 
data for it are available so far. This neglect of triplet correlations has been studied in 



detail for silica melts modeled through the BKS potential 14 1, where c^'^^ was found 
to give noticeable contributions different than in, for example, dense Lennard- Jones 
mixtures. Without triplet contributions, the silica results under atmospheric pressure 
were still qualitatively correct regarding the wave- vector dependence of the correlation 
functions. We expect the same situation for our model potential, with the quality of 
this additional approximation to become better on increasing density. 

From the resulting ^{q,t), an equation similar to Eq. ([T]) allows to calculate 
the a-species tagged-particle density correlations (f)a{q,t) = (exp[i(f- (rait) —raiO))]), 
where r^{t) marks a single tracer particle. The corresponding MCT memory kernel 
readsm^(g,t) - (l/q^) J d^k/{2Trf ^^,^,(g^/9)2c,„,(fc)c„^.(fc)$,,^,(fc, t)</)^(p,i). In 
the limit g ^ 0, the quantity q^m%(q,t) approaches a finite limit which plays the 
role of the memory kernel for the corresponding mean-squared displacement 5r'^(t). 
The self-diffusion coefficients Dq, = limt^oo '^?'Q(i)/(6i) can thus be determined as 
I?„^ = /p dt\m\q-,oq^ml^{q,t). For the numerical solution of Eq. ^ and the 
corresponding equations determining (j)^^(q,t) and fr^(i), we use a wave- vector grid 
with a cutoff of Q = 24/A and grid spacing Ag = 0.1/A. 

MCT describes the slowing down of diffusivity, D^ -^ 0, connected to an 
increase in relaxation times for the ^{q,t), as the coupling described by the Vaa'a" 
increases smoothly through a variation of control parameters (g, T). The divergence of 
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Figure 1. Self-diffusion coefficients of silicon, Dgi (filled diamonds), and oxygen. 
Do (open diamonds), in model silica melts calculated by the mode-coupling 
theory of the glass transition along the indicated isotherms, as functions of density 
Q (top) and as functions of pressure P (bottom). Filled circles: Dgi obtained from 
molecular-dynamics simulation. 



relaxation times defines the MCT critical point Tc{g). No such divergence is observed 
in experiment or simulation, but the scaling laws connected to Tc describing the 
asymptotic shape of the correlation functions and a power-law variation in relaxation 
times are [12]. Hence, the MCT critical point provides a useful, well-defined concept 
to discuss the slow dynamics of glass-forming systems. 
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3. Results 

Figure[T]shows the self-difFusion coefficients Da (a = Si, O) calculated from MCT along 
isotherms as diamond symbols. A maximum is found for all isotherms considered, 
at a density of about go ~ 3.50 g/cm'^ roughly independent of temperature. This 
corresponds to a diffusivity maximum at a pressure of about Pq w 20 GPa, as the 
lower panel of Fig. [1] shows. There, density values have been translated to pressure by 



means of the computer-simulated equation of state [171 lli?|. In both representations, 
the maximum is more pronounced at lower temperatures, indicating that it is a feature 
of slow glassy dynamics. At high temperatures, only the pronounced initial increase 
of Da with pressure remains clearly visible, while the region around the maximum 
becomes rather broad. Note that all the temperatures discussed here are still high 
compared to the conventional glass transition Tg. The maximum in diffusivity also 
corresponds to a minimum in viscosity, or more generally a minimum in structural 
relaxation times. For O diffusion, we find the same behaviour as for Si diffusion, 
while Do is slightly larger than Dsi; the ratio Do/Dsi for T = 3580 K drops from 
roughly 1.78 at g — 2.3g/cm^ to about 1.45 at g — 4.2g/cm^; other isotherms give 
similar behaviour. These ratios are slightly larger than what has been measured in a 
silicate melts with several alkali-oxyde additions [7| at diffusivities of O(10~^^ m^/s). 
Note however that at lower diffusivities, the transport of Si and O will be governed by 
hopping processes with different activation energies, leading to a much larger Do/Dgi. 

The diffusion coefficients from MCT are in qualitative agreement with results 
from MD simulation using the CHIK potential. This is demonstrated for Dsi, 
where Fig. [T] reproduces some simulation data from Ref. [l9| as circle symbols. In 
particular, the density and pressure of maximum diffusivity found in MCT correspond 
well to the computer-simulation results. Already earlier simulation studies based on 
the BKS potential for silica have reported a diffusivity maximum at a density of 
g^^ « 3.5 g/cm^ [9|, llO| indicating that this feature is robust against slight changes 
in the potential. The absolute values of Z?si disagree between MD and MCT. In 
particular at low temperatures and low densities, additional relaxation processes not 
captured in the theory render the divergence of the diffusion coefficient much weaker 
than predicted by the theory, hence the disagreement is most pronounced in this 
regime. Here, computer simulations show a temperature dependence of the Da that 
is well described by Arrhenius laws, which are not reproduced in MCT. 

The strong variation of Da{g) - it increases by orders of magnitude upon 
increasing pressure up to Pq - suggests an explanation in terms of a corresponding 
variation of the critical temperature Tc{g) of mode-coupling theory. The values 
calculated from a bifurcation analysis of the long-time limit of Eq. ([1]) arc shown in 
Fig.[2l Indeed, Tc{g) shows a pronounced minimum around go. For most densities, it 
was possible in the simulation to obtain equilibrated configurations also below this Tc, 
so that the calculation could be based on linear interpolation of a set of state points for 
which S{q) was simulated, and between which the structure factor changes were small. 
Only around go, the critical temperature is too low to equilibrate the system within 
reasonable time scales in MD. Additionally, at g — 4.2 g/cm'^ and temperatures around 
and below T = 2750 K, the simulations showed crystallization, preventing access to 
the liquid regime. Therefore, the Tc values shown for g ~ 3.3g/cm^, 3.5g/cm^, 
and 4.2 g/cm^ are based on extrapolation of S{q) from higher T and have higher 
uncertainty than the rest of the Tc data; the state points at which 'exact' MD input 
for S{q) was used are marked in Fig. [2] by diamond symbols. Translating the (g, Tc) 
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Figure 2. Mode-coupling-theory critical temperature Tc for silica, calculated by 
MCT based on MD-simulated static structure factors with the CHIK potential 
(diamonds), and estimated from the simulation (circles). Open diamonds are 
calculated by extrapolation of MD structure factors S{q) available only for T > T^; 
small triangles indicate state points for which S{q) was available from MD. 



pairs to Tc{P) by the simulated equation of state, the Tc{g) minimum corresponds to 
a similar minimum around Pq (lower panel of Fig. [5]). In particular, the qualitative 
behaviour of Tc{q) and Tc{P) is identical, indicating that thermodynamic features 
governed by the equation of state are not central to understanding the observed 
dynamics. 

Since transport processes slow down dramatically upon approaching Tc, relaxation 
in the silica melt along an isotherm first becomes faster with increasing pressure as the 
distance jT — Tc| increases. Only once the pressure exceeds Pq, this distance decreases 
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Figure 3. Rectification plot for the determination of Tc(^(>) from the MD- 
simulated Si diffusion coefficients Dg;: DJ^'^ with 7 = 2.15 for various densities as 
indicated (symbols). Solid lines are linear fits determining Tc as the intersection 
with the T axis. 



and hence relaxation becomes slower again. For comparison, in the Lennard- Jones 
system, Tci^g) increases monotonically with increasing density [23[ which leads to the 
expected monotonous decrease of diffusivity with increasing pressure. 

To demonstrate the level of agreement between MCT and simulation data, we 
have also estimated T^ from the MD dynamics by fitting asymptotic power laws of 
MCT to the temperature dependence of the diffusion coefficients. The theory predicts 
D ~ |T— Tc|^, and while in principle 7 will depend on density (see below), we anticipate 
that this change will be relatively small and use 7 = 2.15 for all data sets. Figure [3] 
shows exemplary rectification plots from which Tc{q) is determined by linear fits to 
D^l'^ (symbols) in a restricted temperature range not too close and not too far above 
Tci shown as solid lines. Deviations from a linear slope at higher temperatures result 
from preasymptotic corrections, while those at lower temperatures signify non-MCT 
relaxation process ('hopping'). In agreement with the observation made in Fig. [I] the 
latter deviations appear more dominant at lower densities, indicating that the MCT 
description will be better for pressurized silica than under atmospheric conditions. 

Although the T^ determined from simulation, shown as circle symbols in the upper 
panel of Fig. [2l are systematically lower than the Tc calculated within the theory, the 
position of the minimum at qq is in good agreement. There is a slight tendency for 
the MCT result to better agree with the MD-determined T^ at higher densities, in 
line with the expectation that MCT deals quantitatively better with dense liquids. 
Tc at ambient pressure has previously been determined using the BKS potential of 
silica. MD simulations yielded Tc = 3330 K, while MCT calculations without triplet 
correlation contributions gave Tc = 3962 K [1J| . These values are significantly higher 
than our results, although the static structure obtained from the BKS respectively the 
CHIK potential shows relatively small differences [19j . This underlines the importance 
of obtaining static-structure factor input for MCT as accurately as possible. The 
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Figure 4. Nonergodicity parameters F'^{q) calculated within MCT for model 
silica melts, at the critical temperature Tc for densities g = 2.3g/cm'^ 
(solid circles), 3.1g/cm^ (open circles), and 4.2g/cm'' (diamonds). Shown for 
comparison is F'^{q) for a binary hard-sphere mixture (dotted line), see text for 
details. 



inclusion of the static triplet correlation function c^'^^ for the BKS potential led to a 
better agreement in the length-scale dependence of the glass form factors 1^, but 
worsened the agreement for Tc. Whether the same will also hold for the CHIK 
potential, is unclear, but a noticable shift of the Tc values presented in Fig. [2] has 
to be anticipated. 

The variation in Tc shown in Fig.[2]can be understood within MCT by looking at 
specific features of the static structure factors. While at the lower densities, all Sapiq) 
show a pronounced scattering peak at qi « 1.7/ A reflecting Si02 tetrahedra as the 
structural units of the system, this peak vanishes at the expense of the main diffraction 
peak in Sap{q), located at q2 ~ 2.8/A at g = 2.3g/cm'^ and corresponding to typical 
interatomic, rather than inter-tetrahedral distances [l9|. This trend can be identified 
also in the total density-density correlations, S{q) = '^ctS^aX/sSajsiq)- Similarly, it is 
reflected in the glass form factors F{q) = ^^^ limt^oo a;aa;/3$Q/3((7, i), i.e., that part 
of density fluctuations which is frozen in upon crossing the MCT glass transition. 
Figure [4] shows the critical form factors F'^{q) as a function of wave number, evaluated 
at several densities along the MCT transition line Tc{q)- While for gi = 2.3 g/cm'^ (filled 
circles, corresponding to P « 1.25 GPa), two peaks of almost identical height appear 
at the two g- values qi and (72, the first peak has almost disappeared at g — 4.2g/cm'^ 
(diamond symbols, corresponding to P k, 46.2 GPa). At the same time, the peak at 92 
has grown. This evolution reflects a gradual loss of tetrahedral ordering, and a smooth 
crossover to a system showing signatures of a dense liquid. Around the minimum in 
Tc., both contributions prevail, as the F'^{q) for g = S.lg/cm'^ (corresponding to 
P w 13.0 GPa) shown in Fig. [4] demonstrates. Hence we attribute the initial decrease 
of Tc{g) to a loss of chemical short-range order (the tetrahedral structure), and the 
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subsequent increase to an increase in nearest-neighbour cageing. This crossover visible 
in F'^(q) corresponds to a change in mean coordination numbers for the Si atoms in 
Si02 with increasing pressure. MD simulations for example show a gradual crossover 
from predominantly four- fold coordinated Si atoms at g ~ 2.3g/cm'^ to a significant 
number of five- and six- fold coordinated ones [l9| . 

To demonstrate the approach to a frozen structure that resembles one governed 
by packing effects, we show in addition in Fig. 2] the F'^{q) obtained from a binary 
hard-sphere mixture (using the Percus-Yevick approximation for S{q)) with diameters 
chosen as diarge = 1-82 A and dsmaii = 1-46 A and concentration a:iargc ~ 1/3. 
Considering that the covalent radius of Si yields roughly dgj — 2.2 A, the value of diargc 
appears reasonable if one takes into account the complicated interatomic potentials in 
the Si02 melt. In particular, the comparison with the hard-sphere F'^{q) shows that 
besides the growing main peak, the emerging shoulder at g « 5. 7/ A can be attributed 
to excluded- volume effects. It can be anticipated that the F''{q) for the silica melt and 
the hard-sphere mixture further approach each other, as the silica density is further 
increased. Only in the g — » limit the two quantities show no convergence, related to 
the fact that the two systems have rather different isothermal compressibilities. 

A further corroboration for the crossover from a tetrahedral network former to 
an excluded-volume influenced glass comes from the analysis of the MCT exponent 
parameter A. This quantity controls the exponents of the asymptotic expansions near 
the singulary Tc- It is bounded by 1/2 < A < 1, but for common 'fragile' glass formers, 
one usually finds A w 0.75 as, e.g., for the Lennard- Jones model p^, and values 
A > 0.8 for systems where attractive interactions play comparable role to excluded 
volume as, e.g., in the square- well system J24| . The evolution of A with increasing 
density in the present silica model is shown in Fig. [5] Interestingly, it shows values 
larger than 0.8 for the lower densities corresponding to moderate pressures. At high 
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Figure 6. Mean-squared displacement of Si atoms in the modeled silica melt, 
at various densities indicated and along the T = 2750 K isotherm. Dotted lines 
indicate the corresponding quantities at T = Tdg) for g = 2.7g/cm^, S.Sg/cm"^, 
and S.Qg/cm"^. 



densities it systematically decreases towards a value of A « 0.766 sd g = 4.2 g/cm'^. For 
comparison, the value obtained for the binary hard-sphere mixture discussed above, 
A « 0.742, is indicated in Fig. [5] as a horizontal arrow. This suggests that indeed 
the asymptotic dynamic behaviour of the high-density silica melt slowly approaches 
that of a densely packed mixture. Note also that maxima in A, similar to the one 
shown in the Fig. [51 have been argued to arise from an interplay of two different arrest 
mechanisms |24| . consistent with our picture of a gradual crossover with increasing 
pressure. The parameter A in particular determines the exponent 7 for the asymptotic 
divergence of relaxation times or viscosities at Tc] smaller values of A signify larger 7 
and vice versa. 

The non-monotonous variation of Da shown in Fig. [1] leads to a corresponding 
non- monotonous variation in the mean-squared displacement Sr^it) (MSD) upon 
varying density or pressure along an isotherm. We show exemplary results for Sr'^^(t) 
at T = 2750 K and various densities in Fig. [5] The variation at long times, i.e., for large 
displacements Sr^-^ ^ 1 A reflects the change in Da- But also at earlier times, in the 
ps regime, a strong non- monotonous variation in Sr^it) remains. Since asymptotically 
close to Tc, the plateau visible in the MSD is a measure of the localization length of 
the individual particle, it is tempting to read off from the intermediate-time window 
in Fig. [6] a Si-localization length that shows an apparent change of almost a factor 2 
as a function of density. This is what is also observed in the MD simulation [19| . In 
fact, the value of, say, 5rgj(lps) first increases with increasing density or pressure: Si 
particles are on that time scale less localized for higher pressure, reflecting a change 
in the microscopic dynamics that becomes less dominated by the strong localization 
within the tetrahedral network. However, this is not the localization that is responsible 
for the MCT glass transition, and which can only be read off from the MSD sufficiently 
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close to Tc- To this end, we show as dotted hnes in Fig. [S]the Sr^^{t) a.t T — Tc{g) 
and observe that they agree closely over the entire density range presented, indicating 
that the cage localization length rgj w 0.15 A changes only weakly as a function of 
pressure in this system. If observed far from Tc, large preasymptotic corrections that 
are particularly dominant for the MSD [25| give rise to the apparent shift in r'^. 

For the O atoms, the same qualitative trend as discussed in connection with Si 
holds, although the variation of Tq is slightly larger. The ratio of localization lengths 
6 = ?'o/'"si monotonically decreases from S « 1.25 a.t g = 2.3g/cm'^ to about S w 1.05 
at g = 4.2g/cm^, indicating that the role played by the Si and the O atoms in the 
dynamics assimilates at large pressure. Overall, both localization lengths agree with 
Lindemann's criterion for melting, stating that r^ « O.Ida. 

4. Conclusions 

We have demonstrated that MCT, together with computer-simulation input for the 
equilibrium liquid structure, reproduces a peculiar change in the dynamics of a 
pressurized silica melt: upon increasing pressure, atomic-scale transport as monitored 
through, e.g, self-diffusion coefficients, first becomes faster. At a pressure around 
20 GPa, a maximum in diffusivity occurs, and at still higher pressures, transport 
starts to slow down with increasing pressure. This is in broad agreement with 
previous simulation data on model silica melts and with experiments on various silicate 
mixtures. 

MCT explains the diffusivity maximum in silica melts as arising from a gradual 
change in the static structure on mesoscopic length scales, where contributions 
connected to tetrahedral ordering at wave number q « 1.7/ A become continuously 
less pronounced at the expense of contributions on the length scale of the Si-Si atom 
nearest-neighbour distance, q « 2.8/A. This is a feature found only in silica and similar 
network-forming melts, whereas in simpler dense liquids such as the Lennard- Jones 
liquid, the nearest-neighbour contribution in S{q) remains dominant at all densities 
or pressures. Hence, the latter 'fragile' liquids do not show diffusivity maxima. As an 
additional difference, the variation of Tc{P) with pressure in Lennard- Jones liquids is 
dominated by thermodynamic contributions (viz., a strong variation in the equation 
of state) arising from the presence of a gas-liquid spinodal [2j| , while in silica melts 
it is inherent to the slow glassy dynamics and not governed by the equation of state. 

The qualitative correctness of MCT predictions for pressurized silica melts is 
remarkable, since most experimental observations concern temperatures well below 
the MCT Tc, where the theory in its present form is not applicable. Nevertheless, the 
qualitative physical mechanisms responsible for the anomalous pressure dependence 
of transport coefficients in this strong glass former seem to be there already at far 
higher temperatures, and captured in the MCT approximation. A similar conclusion 
might apply to the distinction between 'strong' and 'fragile' behaviour of the viscosity 
r]{T) around Tg. While MCT cannot be applied there, it does yield for siHca a 
gradual crossover of the exponent 7 governing the initial increase in viscosity above 
Tc, changing from a slower increase at ambient pressure to a stronger one (more akin 
to a 'fragile' liquid) at high pressures. 

Based on this observation, we suggest that MCT can be used to investigate in 
more detail the connection between 'strong' glass formers with network-like structures 
at low pressures and 'fragile' ones which are characterized by a dense arrangement of 
constituent atoms. 
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